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Abstract 

The cross section for the photoproduction of large-pT <7/b rnesons at HERA is 
calculated at next-to-leading order, adopting a perturbative approach to describe 
the fragmentation of charm quarks and gluons into J/p mesons. We treat the 
charm quark according to the massless factorization scheme, where it is assumed 
to be one of the active flavours inside the proton and the resolved photon. We 
present inclusive distributions in transverse momentum and rapidity, including the 
contributions due to direct and resolved photons. The importance of the colour- 
octet components of the J/ip wave function, which contribute to the fragmentation 
process, is emphasized. In addition to prompt J/ip production, we consider also the 
production of Xcj states followed by radiative decays to J/ijj mesons, both in the 
colour-singlet and colour-octet channels. 
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1 Introduction 


The production of heavy quarkonia in high-energy collision has conventionally been de¬ 
scribed within the colour-singlet model (CSM) |]^. In this model, a non-relativistic ap¬ 
proximation is used to describe the binding of the heavy-quark pair, produced via parton- 
fusion processes, into a quarkonium state. The heavy-quark pair is projected onto a 
physical quarkonium state using a colour-singlet projection and an appropriate spin pro¬ 
jection. Experimentally, J/il) production is measured most easily through the detection of 
the leptonic decays. Then, the inclusive J/tf) yield consists of the prompt J/'0 production 
and the production of the cc ^Pj states Xcj {J = 0,1, 2), which decay electromagnetically 
into the lower-lying J/'ijj state. The latter must also be taken into account if one wishes 
to compare theoretical predictions with experiments where Xcj production could not be 
measured separately. The CSM has been successfully applied 0 to explain large-pr J/4’ 
production at the relatively low energies (a/s = 63 GeV) of the CERN ISR experiments. 
After the advent of the CERN SppS collider with y/s = 630 GeV, it was found 0 that 
also the production of bottom quarks and their subsequent weak decay into J/pj mesons 
is an important source of J/pj production. In the UAl experiment, it was not possible to 
separate the bottom-quark contribution. 

Later, the CDF and DO collaborations, when measuring J/pj production at the Fermi- 
lab Tevatron (a/s = 1.8 TeV), were able to subtract the bottom-quark contribution from 
the total yield. Since also the yield of the Xcj mesons decaying into J/pj mesons could 
be separated, a measurement of prompt J/p) production at large pt was possible. (We 
shall use the term prompt instead of direct to denote all J/"/ mesons not originating from 
decays, i.e., R-meson or Xcj decays.) The cross section of prompt inclusive J/'/ produc¬ 
tion measured by CDF and DO turned out to be an order of magnitude larger than 
the prediction of the CSM. It is well known that, at large pt, the dominant mechanism 
of light-quarkonium production is via the fragmentation of gluons and quarks, while the 
parton-fusion process is suppressed by powers /p^, where m is the mass of the quarks 
bound in the quarkonium state. So, the fusion process, which is a higher-twist contribu¬ 
tion, can contribute signihcantly only at extremely small pt, which for the light quarkonia 
is in the soft region outside the perturbative regime. For J/pj production, this region, due 
to the larger charm-quark mass rUc, extends to larger pt- But eventually, at px S> me, 
the fragmentation process, although of higher order in as compared to the fusion process, 
can exceed the fusion process. In fact, it was pointed out by Braaten and Yuan [Q that, 
in addition to the parton-fusion contributions described by the CSM, fragmentation of 
gluons and charm quarks is an important source for large-p^ mesons in hadron-hadron 
collisions. 

J/ifj production in high-energy ep collisions at HERA is dominated by photoproduc¬ 
tion, where the electron is scattered by a small angle producing photons of almost zero 
virtuality. First experimental results in the kinematical region 30 GeV< W < 150 GeV, 
where W is the yp centre-of-mass (CM) energy, have been presented by the HI collabora¬ 
tion 1^. In this work, the inelastic process is analysed for the hrst time. Most of the events 
come from the elastic and diffractive processes. They can be separated by measuring the 
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J/ip energy spectrum in terms of the inelasticity variable z = Pp -pj/p/pp with pp, 
and pj/^ being the four-momenta of the proton, the photon, and the J/ifj meson, respec¬ 
tively. For elastic/diffractive events, one has z ^ 1. A clean sample of inelastic events can 
be obtained by selecting events with 5 < 0.9. Due to the limited integrated luminosity in 
the 1994 HERA operation, the inelastic J/'ip production cross section da/dp^ could only 
be measured for pf, < 6 GeV^ 0. In this range of p^, the inelastic production proceeds 
predominantly through the fusion mechanism in the GSM. In fact, the experimental pt 
and 5 distributions are very well described, as for both shape and absolute normalization, 
by the GSM next-to-leading-order (NLO) calculation of Kramer 0, using up-to-date NLO 
parton density functions (PDF’s) of the proton. 

As is well known, high-energy photoproduction proceeds through two distinct mech¬ 
anisms: direct photoproduction, where the photon couples directly to the quarks, and 
resolved photoproduction, where the photon interacts through one of its partonic compo¬ 
nents, quarks or gluons, via hard scattering with the partonic constituents of the proton. 
For J/'ll) photoproduction within the GSM, the resolved-photon process is 

1 + 9-1 + 9 p^ cc[l,^Si]+g ^ J/ij + g. (1) 

This means that the gluon coming from the proton interacts with the gluon from the 
proton to create a ^Si colour-singlet cc pair, which subsequently fuses into a meson. 
This process has been analysed with the result that its contribution to the overall cross 
section in the GSM is small everywhere, except in the very-low-5 region. Gontributions 
from the production and radiative decay of Xcj states are of similar magnitude; see Ref. 
for a review of these and other resolved-photoproduction mechanisms within the GSM. 

Within the factorization approach, J/^ production can also proceed through colour- 
octet cc pairs. They appear for direct photoproduction as well as for resolved photopro¬ 
duction. Golour-octet conhgurations which contribute to inelastic direct J/ijj photopro¬ 
duction are produced through the following subprocesses: 

X + g ^ cc[8,%] +g, 

X + 9 ^ cc[ 8 ,^Si] +g, 

X + g ^ cc[ 8 ,^Pj]+ g, (2) 


where J = 0,1,2. Other contributions, such as light-quark initiated ones, are strongly 
suppressed at HERA energies and can be safely neglected. Adopting the colour-octet 
matrix elements obtained through hts to prompt J/ijj data from the Tevatron §], one hnds 
that colour-singlet and colour-octet contributions to inelastic direct photoproduction are 
of comparable size fl^ • In the case of resolved photoproduction, the following colour-octet 


processes contribute: 


9-1 + 9p^ cc[8, ^S'o] + g, 

9 - 1 + 9 p^ cc[8,^S'i] + g, 

9 -i + 9 p^cc[ 8 ,^Pj]+g. (3) 
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These are the same processes that contribute significantly io J/ip production in pp colli¬ 
sions at the Tevatron j^]. Their contribution to inelastic J/ip photoproduction has also 


been calculated []^ and found to be non-negligible. In fact, it is comparable to the cross 
section of colour-singlet direct photoproduction 0 - The distributions in the inelasticity 


variable 5 are, however, quite different. The resolved-photon cross section (colour-singlet 
plus colour-octet) is dominant at small z, whereas the direct-photon contribution is large 
towards z = 1. 

For Xcj production, the colour-octet channels are much more important. In the case 
of direct photoproduction, the production of colour-singlet Xcj states via 


7 + g ^ccli, ¥j] + g ^ Xcj + g 

is forbidden to leading order in so that the colour-octet processes 


(4) 


7 + g 

7 + q(q) 


cc[ 8 ,^S'i] +g, 
cc[ 8 ,^^i] +g(g), 


(5) 


where the colour-octet ^Si cc pair fuses into a physical Xcj particle, are the dominant 
ones. Nevertheless, using the matrix elements obtained through fits to the Tevatron data 
PJT^JI3| , the contribution of the channels is still small compared to the dominant CSM 


contribution to J/ip production [|^]. The cross sections for resolved photoproduction of 


Xcj mesons in the colour-octet channels have not yet been estimated. 

All these results refer to J/ip production at pt values of order rric, where we may expect 
the higher-twist contribution to be dominant. At sufficiently large pt, the fragmentation 
process is supposed to be the dominant mechanism oi J /ip photoproduction. First esti¬ 
mates of the fragmentation contribution to J/ip photoproduction have been presented in 
Ref. 1^; see also Ref. for earlier work within the CSM. These authors considered, 
at leading order (LO), the fragmentation of charm quarks and gluons produced through 
the direct-photon processes 7 + g ^ g + g and 7 -|- —*• c - 1 - c. They did not take into 

account resolved photoproduction, which is certainly important, too. In the case oi J /ip 
production in pp collisions at the Tevatron, it is the analogous mechanism which makes 
the dominant contribution 
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It is the purpose of this work to study J/ip photoproduction at large pt including 
all direct- and resolved-photoproduction channels at NLO. We shall focus our attention 
on the production via fragmentation. We shall consider the fusion mechanism only for 
comparison; for simplicity, we shall only include its dominant LO contribution, which 
arises from direct photoproduction in the CSM. The superposition of the NLO cross 
sections for the fusion and fragmentation mechanisms is left to future work. 

The framework for our calculation of the NLO fragmentation cross section is similar 
to the one employed recently for the calculation of inclusive D*^ photoproduction . 
At large pt, one may neglect the mass of the charm quark. The collinear singularities 
corresponding to the asln(p|./m^) terms of the massive-charm scheme, which is used, 
for example, for predictions at pr^^ci are then absorbed into the charm-quark PDF’s 
and the fragmentation functions (FF’s) in the same way as for the lighter m, d, and s 
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quarks. The FF’s for the fragmentation of charm quarks and gluons into J/'^ mesons 
must be known for this approach at the scale of production, which is usually taken to 
be /i = rriT, where rriT = is the transverse mass of the J/iIj meson with 

mass Mj/^ ^ 2mc- These FF’s can be computed at the starting scale po = Mj/^ in 
nonrelativistic QCD adopting as nonperturbative input the colour-singlet and colour- 
octet transition-matrix elements, which also enter the fusion approach. At larger scales 
fi S> Mj/^, the large logarithms in are resummed using the usual Altarelli-Parisi 

evolution equations. At the starting scale, only the fragmentation of gluons and charm 
quarks need to be included, while the light-quark contributions are negligible. The latter 
will be generated at higher scales via evolution. 

The outline of our work is as follows. In Section 2, we shall shortly describe the 
formalism for the perturbative FF’s and discuss the transition from massless to massive 
factorization used also in our previous work |^. In Section 3, we shall present our 
NLO predictions for the various fragmentation contributions to the cross section of J/ip 
photoproduction in ep collisions at HERA as functions of pt and rapidity y, and compare 
them with the LO CSM results based on the cc fusion mechanism. Specihcally, we shall 
distinguish between direct-photon and resolved-photon, prompt and non-prompt (only via 
Xcj states), colour-singlet and colour-octet fragmentation production. Our conclusions 
will be summarized in Section 4. 


2 Charmonium fragmentation functions 


In this section, we shall describe the assumptions underlying the massless approach, which 

For the reader’s conve- 
As already mentioned 


is very similar to that in our work on inclusive D* production 0 


nience, we shall repeat here the basic ideas developed in Ref. [|I^ 


in the Introduction, two mechanism contribute to the photoproduction of charm quarks 
in ep collisions: {i) In the direct photoproduction mechanism, the photon couples directly 
to the quarks, which besides the massless u, d, and s quarks, also include the massless 
charm quark. In this case, at least in the LO processes, no spectator particles travel along 
the momentum direction of the photon, [ii) In the resolved photoproduction mechanism, 
the photon splits up into fluxes of u, d, s, c quarks and gluons, which then interact with 
the partons coming from the proton leading to the production of quarks, including charm 
quarks, and gluons at large pt- The contributing parton-level processes are the same as 
in the case oi J /'ip production in hadron-hadron collisions. The quarks and gluons coming 
out of the photon are accompanied by a spectator jet, which also travels in the photon 
direction. Therefore, the yp cross section depends not only on the PDF’s of the proton 
but also on those of the photon. The main difference with respect to the usually consid¬ 
ered massive-charm scheme, which to our knowledge has not yet been worked out for J/ip 
production via fragmentation, is that the charm quark also contributes via the PDF’s of 
the proton and the photon, he., charm is already an active flavour in the initial state. 
Compared to production, where the mesons are almost exclusively produced via 
charm-quark fragmentation, m. J /ip production the fragmentation process g — >■ J / ip + X is 
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very important as well, as we shall see in the next section. The massless-charm approach 
is justihed if a large scale is governing the production process. In our case, this is the 
transverse momentum of the J/'ip meson, with pt 3> Mj/^p. In this region, non-singular 
mass terms are suppressed in the cross section by powers of ijIc/pt- The important mass 
terms appear if the virtuality of the charm quark is small. This occurs in the initial state 
if the charm quark is emitted from the proton or photon, and in the hnal state if the 
partons emit a cc pair. These two contributions lead to the ln(py/m^) terms in the 
massive scheme. In our massless approach, they are summed into the scale-dependent 
proton and photon PDF’s and J/ip FF’s. 

To obtain the hnal J/ip production cross sections in NLO, the following steps are 
taken, (i) The NLO hard-scattering cross sections for the direct- and resolved-photon 
processes are calculated in the massless approximation with nj = A active Favours |T7[ . 
The collinear singularities are subtracted according to the MS scheme. Since the charm 
quark is taken to be massless, the singularities from its splittings are subtracted as well. 
{ii) The charm quark is accommodated in the PDF’s of the proton and photon as a light 
Favour. The Fnite mass of the charm quark is taken into account by including it in the 
evolution of the PDF’s in such a way that these are only non-vanishing above a scale 
set by its mass. {Hi) The FF’s characterize the hadronization of the massless partons, 
including the charm quark, into cc bound states, he., the J/ijj and Xcj mesons. Similarly 
to the fragmentation into light mesons, these FF’s are non-perturbative input and must 
be determined by experiment, for example from the cross section of inclusive J/V’ and XcJ 
production in hadron-hadron collisions. This information determines the FF’s at some 
starting scale, /iq. 

An alternative methods is to calculate, within perturbative QCD, universal starting 
conditions for the FF’s at a scale po of order This allows one to include in a 

straightforward manner the colour-octet contributions with parameters related to the 
ones obtained from the analysis of colour-octet fusion. 

The structure of the perturbative FF’s is based on the general factorization analysis of 
the production and decay of heavy quarkonia by Bodwin, Braaten, and Lepage |^. This 


factorization formalism allows for the FF’s of charmonium production to be factorized 
into short-distance coefFcients that describe the production rate of a cc pair within a 
region of size 1/mc, and long-distance factors that contain the non-perturbative dynamics 
responsible for the formation of the bound state H containing the cc pair. The FF, 
at factorization scale p, for a parton a to fragment into the charmonium state H with 
longitudinal momentum fraction z can be written as 




( 6 ) 


where n labels diFerent colour components (singlet or octet) and spin quantum numbers 
to be speciFed later. are the local four-fermion operators deFned in Ref. 0 in 
terms of the Felds of nonrelativistic QCD. The short-distance coefFcients d'2{z,fi) can be 
computed at the starting scale po using perturbation theory in as{po). The dependence 
on the observed charmonium state H only appears in the long-distance factors |0). 
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As input for our analysis, we need the short-distance factors d'^{z, /xq) and matrix elements 
(0|C>,f |0) for the following transitions: g —> J/x/', g —> Xcj, c —>• J/V’, and c —> Xcj, which 
we shall take from the literature. In the case of the gluon FF’s into J/^jj and XcJ mesons, 
we use the results of Refs. and |^, respectively. The corresponding charm-quark FF’s 
may be found in Refs. and |^, respectively. 


In order to understand our hnal results, it is necessary to know the relative importance 
of the various terms labelled n in Eq. (P). On the one hand, this is determined by the 
magnitude of the short-distance factors and, on the other hand, by the relative 

size of the various matrix elements (0|(9^|0). The size of the latter may be estimated by 
how they scale with the relative velocity v of the charm quarks inside the charmonium 
state. The magnitude of the d'^{z^ jj) functions is controlled by the lowest order of that 
appears in d'^{z,fi). Thus, to determine the relative importance of the terms in Eq. (H), 
one needs to know the scaling in v of the matrix elements and the order in of their 
coefficients. 

Let us start with the gluon FF into J/V’, which is a ^Si cc state. Keeping only the 
leading terms, Eq. takes the form 


m-i 


m-i 


where d[-’^^^Hz) and d[-’^^^ 


are dimensionless functions listed in the Appendix. 


The hrst term in Eq. ( 0 ) corresponds to the colour-singlet contribution leading in v. Its 
matrix element (0|(9/’^’^(^S'i)|0) is proportional to the probability for the formation of a 
J/il) meson from a pointlike cc pair in a color-singlet ^Si state. It is related to the J/'0 
radial wave function at the origin, Rjp(O), by |2^ 


{0|C>i'«(=>Si)|0) = 


( 8 ) 


up to corrections of relative order and may be extracted from the partial width of the 


leptonic decay J/'ip ^ . Following Ref. we use (^51)10) = 1.13 GeV' 


which includes the 0{as) correction to the leptonic width. The leading contribution to 
the short-distance coefficient comes from the parton process g* —>■ ccg, where 

g* denotes a virtual gluon, and is of order af. The matrix element (0|(!l/^’^(^S'i)|0) is of 
order v^, so that the contribution to the gluon FF is of order All other matrix 

elements in the expansion in Eq. (0) are of higher order in and would, therefore, be 
neglected in the colour-singlet approximation. However, in the charmonium system, 
is not actually small; typically one has ~ 1/4, so that not all matrix elements that 

are formally suppressed by powers of may be neglected in practice. Of particular 

importance is the matrix element (0|(!lg^^(^S'i)|0) because it has a short-distance factor of 
order which which arises from the parton process g* —> cc. The corresponding matrix 
element (0|C>8^’^(^S'i)|0) is of order We take its numerical value to be (0|C>8'^’^(^S'i)|0) = 
0.014 GeV^ [2^. Thus, the suppression factor is compensated by the factor 1/a^. 
In addition, the different ^ dependence of d\^’^^^\z) as compared to the colour-singlet 
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short-distance coefficient enhances the relative importance of the colonr-octet 

contribntion. 

The glnon FF’s into the ^Pj charmonia Xcj have been worked ont in Ref. [R] and take 
the form 

= ^%i4i*P.l(j)(0|c>w(3p_,)|0) + ^£Mdp.l(..)(0|O*-3,35^)10), ,9) 


m'i 


m-i 


The colonr-singlet matrix element (0|Clf‘'‘^(^Pj)|0) is related to the derivative of the non- 
relativistic radial wave fnnction of the P-wave states at the origin, P^^(O), by 
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(o|c>f-^(¥j)|o) = (2J + 1)—|p; 


27r' 


Xc 


( 10 ) 


np to corrections of relative order This parameter can be determined phenomenolog¬ 
ically from the annihilation rates of the Xcj mesons. Since both matrix elements are of 
the same order in v and the colonr-octet term has one power of ag less, it is expected to 
dominate. According to Refs. and |^, the matrix elements take the valnes 

1 


(0|Of°(^Po)|0) = 


(0|C>f‘^('^Pj)|0) = 0.0862 GeV^ 


2J+ 1 

(0|0|=°(35i)|0) = ^^(0|O8^‘=‘^(3Ri)| 0) = 0.0076 GeV^ 


( 11 ) 

( 12 ) 


In the evalnation of Eq. (|^), we have used 2mc = Mj/^ = 3.09688 GeV |^. The value of 


the colour-octet matrix element in Eq. (R) is larger than the most recent value obtained 
from P-meson decays 0. 

We now turn to charm-quark fragmentation. The LO formula for the c ^ J/ip FF 
has been found in Ref. 0 


and takes the form 


m-i m-i 


where and may be found in the Appendix. While, in the case of 

g the v'^ suppression of (0|C>8^^(^S'i)|0) relative to (0|C>/^^(^S'i)|0) is compensated 

by the fact that the colour-octet coefficient is enhanced by a factor l/a^ relative to the 
colour-singlet one, such a compensation does not occur in Eq. (^). In fact, the second 
term in Eq. (0) only amounts to about 0.12% of the hrst one. 

The c —>• Xcj FF’s are given by [^] 


= £lW4nf3i(.)(0|of-(3p^)|0) + 5lW4aA.i(.)(0|o3-(3s^)l0)_ , 14 ) 

m° m-i 


where {z) and {z) are specihed in the Appendix. As pointed out in Ref. 

the colour-octet term in Eq. (|0) is minuscule. Although the minimum invariant 
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mass 









































of the fragmenting charm quark is SrUc, for consistency with the gluon FF’s, we choose 
/xo = 2mc in Eqs. (|T^ and (JT^). 

For the calculation of the J/'ip production cross sections, we need the FF’s at the 
factorization scale fi S> Mj/^. Then, large logarithms in appear, which have to 

be resummed. This is achieved by using the Altarelli-Parisi equations. 







(15) 


where 


HP {z,as{^)) = 


(/^) 

27r 


A°’^T) + 


2tt 




(IT), 


ba 


z) + 0 {a 


(16) 


are the timelike splitting functions of parton a into parton b. In our NLO analysis, 
we include the splitting functions through order [^, while in our LO analysis we 
truncate Eq. (|T^) after the first term. The initial conditions D^{z,^o) for a = g,c and 
H = J/'ip,Xcj have been specified above. Due to charge-conjugation invariance, we have 
D^{z,fio) = D^{z,fio). The initial light-quark FF’s are set equal to zero. They are 
generated at larger scales /i via Eq. ( 0 )- However, we expect their effect to be small. 
The splitting functions pI;P {z,as{iJ.)) are for massless quarks throughout. 

In Ref. an improved version of Eq. (pISD , which respects the phase-space constraint 
Daiz,g,) = 0 for ^ < M^//x2, 


has been suggested. In Ref. |26 


it was found that the 

numerical effect of this improvement is not crucial at large scales. Since we are mainly 
interested in the fragmentation production oi J/ip mesons at large pt, for simplicity, we 
shall stick to Eq. ([T5|) for the time being, leaving the implementation of the generalized 
evolution equations to future work. 

As in our recent work on D*^ photoproduction [jT^, we start from the NLO hard- 
scattering cross sections calculated in the MS scheme with massless flavours, and slightly 
modify the factorization scheme for the collinear singularities associated with hnal-state 
charm quarks. Specihcally, whenever final-state collinear singularities are subtracted from 
the hard-scattering cross sections, we substitute P^^’R{z) \n{s/jjp) P^^’P{z) ln(s/p^) — 

dca{z), where P^^’P{z) are the LO timelike a ^ c splitting functions of Eq. (pTl) and 
the dca{z) functions may be found in Ref. [|T^. Here, a = g,c,q, where q stands for 
the hrst three quark flavours. With this procedure, the factorization of the hnal-state 
collinear singularities associated with the charm quark is adjusted so as to match the 
hnite-ruc calculation. This is equivalent to the matching approach proposed in Ref. 
between the massless-charm calculation in connection with the perturbative FF’s and the 
massive-charm calculation without FF’s. 


3 Results 

This section consists of three parts. Firstly, we shall specify our assumptions concerning 
the proton and photon PDF’s as well as the equivalent photon approximation. Secondly, 
we shall present our predictions for the particular range of yp CM energies W that is used 
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in the HI analysis of inclusive J/ijj production P, namely, 30 GeV <W< 150 GeV. In 
this subsection, we shall compare the results for the various channels of interest: prompt 
J/ijj photoproduction versus non-prompt J/'ip production originating from Xcj photopro¬ 
duction with subsequent XcJ — J/ip + 'y decay; colour-singlet channel versus colour-octet 
channel; direct photoproduction versus resolved photoproduction. In the third part, we 
consider the high-lH range, 150 GeV < W < 280 GeV, where the fragmentation contri¬ 
bution is greatly increased relative to the fusion contribution. 


3.1 Input information 

For the calculation of the cross section dPa/dyiabdpT, we adopt the present HERA condi¬ 
tions, where Ep = 820 GeV protons collide with = 27.5 GeV positrons in the laboratory 
frame, so that -y/s = EpE^, = 300 GeV is available in the GM frame. We take the ra¬ 
pidity yiab to be positive in the proton flight direction. The quasi-real-photon spectrum 
is described in the Weizsacker-Williams approximation by the formula 


= 


a 


1 + (1 -a;)^ 


n2 
In 


X 


Q 


2 

min 


+ 2m^x 




2 

max 



(17) 


where x = E^/E^ = W^/s, Ql^iin — '^ea^^/(l — x), a is Sommerfeld’s hne-structure 
constant, and rrif. is the electron mass. We consider the case where the hnal-state electron 
is not tagged, so that ^ GeV^. The yp energy interval 30 GeV < IV < 150 GeV 

(150 GeV < IV < 280 GeV) considered in Section 3.2 (3.3) corresponds to 0.010 < x < 
0.249 (0.249 < x < 0.869). 


We work at NLO in the MS scheme with Uf = A flavours. As for the proton and 
photon PDF’s, we use set GTEQ4M [^, with A^ = 296 MeV, and set GRV HO 
after converting it from the DIS.y scheme to the MS scheme. We identify the factorization 
scales associated with the proton, photon and hnal-state hadron and collectively denote 
them hj Mf. For Mf and the renormalization scale p, we choose y = Mf = rriT, where 
rriT is the J/ip transverse mass dehned above. We take the starting scale of the FF’s 
to be po = 2mc. We calculate as{y) from the two-loop formula with A^ equal to the 
value used in the proton PDF’s. Unfortunately, the NLO correction to the perturbative 
coefficient functions discussed above, which, strictly speaking, ought to be included as 
well, are not yet available. However, in analogy to the genuine NLO corrections io J/ip 
photoproduction at HERA in the GSM, which are only at the 20% level, these are likely 
to be moderate as well. 

The LO results which enter the QGD-correction {K) factors are calculated consistently, 
he., using the one-loop formula for as(/i), the GTEQ4L proton PDF’s, with A^ = 
236 MeV, and the GRV LO photon PDF’s, and evolving the FF’s with the LO splitting 
functions. In all calculations, we employ the numerical values of the matrix elements 
specihed above. 
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3.2 Results for the low-VF range 

In this subsection, we investigate the NLO cross section (fa/dyiabdpT of inclusive J/'0 
photoproduction at HERA, integrated over 30 GeV < W < 150 GeV as in the HI 
experiment 0 . In order to analyze the pt dependence of the cross section in the central 
region of the detector, we integrate d^ajdyiahdpT over —1.5 < yiab < 1- The result 
is displayed in Figs. |^(a)-(c) for direct and resolved photoproduction and their sum, 
respectively. In each case, we show how the total fragmentation contributions (solid lines) 
are composed of the respective prompt colour-singlet (dashed lines), prompt colour-octet 
(dot-dashed lines), and non-prompt (dot-dot-dashed lines) contributions. For comparison, 
we also show the LO direct-photon contribution of the GSM (dotted lines). We have to 
bear in mind that the latter is enhanced by typically 20% due to genuine NLO QGD 
corrections to the parton-fusion cross section [0. However, we implicitly include the bulk 
of the QGD corrections, namely, those to the leptonic J/'0 decay width, which affect 
the value of (0|(!l/^’^(^S'i)|0) extracted from experiment [^. For consistency with the 
fragmentation contribution, we also choose p = Mf = rriT in the fusion contribution. 
We do not take into account the resolved-photon GSM contribution, which appreciably 
contributes only at very low z, i.e., at low pt and large yiab H]- A dedicated study of the 
low-5 range, which we intend to do in the future, should also consider this contribution. 

As anticipated in the Introduction, both in the direct- and resolved-photon channels, 
the fragmentation cross sections fall off less steeply with p^ than the fusion cross section. 
In all cases, we observe that the prompt colour-octet contribution is most significant and 
dominates the large-p^ behaviour. The non-prompt cross section is always smaller than 
the prompt cross section. In the sum of the direct- and resolved-photon channels, their 
ratio ranges from 2.8 at = 3 GeV to 4.3 at pt = 15 GeV. We emphasize that, in 
the low-pT range, the resolved-photon component, which was previously neglected [|T^, is 
considerably more important than the direct-photon component, by a factor of 18 (2.3) 
at Pt = 3 GeV (10 GeV). The two components cross over at pt ~ 15 GeV. Notice also 
that, even for Pt^Q GeV, the resolved-photon component is comparable to the fusion 
contribution; for p^ > 7 GeV, it exceeds the fusion contribution. In the resolved-photon 
contribution, the prompt colour-singlet part is negligibly small compared to the prompt 
colour-octet part. This is caused by the dominance oi g J/'iIj+X in the resolved-photon 
cross section and the suppression of colour-singlet gluon fragmentation mentioned above. 
In the direct-photon cross section, this hierarchy is reduced since c—*• J/'iIj + X is more 
important. 

In general, rapidity distributions allow for a more specihc investigation of the indi¬ 
vidual contributions than transverse-momentum distributions. In Figs, ^(a)-(c), the yiab 
dependence of d'^a/dyiabdpT is displayed at p^ = 5 and 10 GeV for the same cases as in 
Figs. 2(a)~(c). While the direct-photon component peaks in the same yiab region as the fu¬ 
sion contribution, the resolved-photon cross section dominantly contributes in the forward 
direction, with a maximum at about yiab = 2.5. As might be inferred from Figs. |^(a)-(c), 
the prompt colour-octet contribution is dominant; in fact, this is true over the full yiab 
range. We observe that the yiab spectra of the prompt colour-singlet channel, although 
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suppressed, have a shape similar to the fusion contribution and are peaked at signihcantly 
smaller yiab values than the other fragmentation channels. At large pt, the dominant J/V’ 
production mechanism in pp collisions at the Tevatron is colour-octet gluon fragmentation 
|30| . In order for HERA experiments to probe this component, it will be necessary to 


take data in the very forward direction, where the resolved-photon colour-octet contri¬ 
bution is most signihcant. It should be clear that the colour-octet contributions shown 
in Figs. §(a)-(c) are proportional to the corresponding matrix elements specihed in Sec¬ 
tion 2, which are extracted from hts to Tevatron data. So, in order to substantiate (or 
question) the huge prompt colour-octet J/V’ contribution that is apparently indispensable 
to reconcile the Tevatron data with theory, it would be desirable to collect data in the 
very forward region at HERA, which has not yet been explored by the HI and ZEUS 
collaborations. 


As mentioned in the Introduction, the present analysis extends the previous work [|14 
on J/pj fragmentation production at HERA in two important respects. For one thing, 
we include the resolved-photon contribution, which was not mentioned in Ref. |]T^. On 
the other hand, we work at NLO, while the analysis of Ref. proceeds at LO. Above, 


we have stressed the signihcance of the resolved-photon contribution, which is nicely il¬ 
lustrated in Figs. |l](b) and @(b). In the following, we shall assess the impact of the NLO 
corrections on the various fragmentation contributions. To this end, we consistently re¬ 
peat the analysis of Figs. |l|(a)-(c) at LO, as described in Section 2, and plot the resulting 
NLO to LO ratios in Figs, ^(a)-(c). In the case of direct photoproduction considered in 
Fig. |^(a), the prompt colour-singlet contribution has an approximately constant K factor 
of about 0.6, while the prompt colour-octet and non-prompt contributions have dimin¬ 
ishing K factors only at pt values below 6.5 and 4, respectively. In the case of resolved 
photoproduction addressed in Fig. |^(b), the K factors of the various channels are always 
signihcantly larger than unity, and the overall K factor ranges between 2 and 2.4. In both 
the direct- and resolved-photon cases, the K factors of the non-prompt contributions take 
extraordinarily large values in the upper pt range. Detailed investigation reveals that 
this enhancement essentially originates from the lower edge of the kinematically allowed 
Hiah range. However, since, at large pt, these contributions are greatly suppressed relative 
to the prompt colour-octet contributions, this is inconsequential for the superposition of 
all channels, which is experimentally observable. In fact, the overall K factor, which is 
shown as the solid line in Fig. ^(c), is well behaved, ranging between 1.6 and 1.8 in the 
Pt range considered. 

At this point, we should comment on the LO analysis of direct J/pj photoproduction 
via fragmentation at HERA reported in Ref. fl^. Similarly to the present study, the 
authors of Ref. 0 included the prompt, non-prompt, colour-singlet, and colour-octet 
contributions. However, they did not convolute the yp cross section with the Weizsacker- 
Williams distribution, but considered hxed photon energies instead. There is an overall 
factor of 18 missing on the right-hand side of Eq. (3) in Ref. [M|, which is to describe the 
fusion cross section. In the case of fragmentation production, a meaningful quantitative 
comparison is impossible, since the input values of the poorly known colour-octet matrix 
elements are not to be found in Ref. [ITl . 
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One should bear in mind that all results presented in this subsection refer to the low- 
W range 30 GeV < W < 150 GeV. We expect that the resolved-photon contribution will 
be enhanced relative to the fusion contribution if W is increased. In the next subsection, 
this issue will be addressed in some detail. 


3.3 Results for the high-VF range 


In order to test the predictions for J/'ip photoproduction via fragmentation at HERA that 
follow from the theoretical analysis of the Tevatron data, it is useful to hnd regions of 
phase space where the fusion mechanism is suppressed. As we have seen in Figs. @(a)- 
(c), this happens, for instance, in the very forward region of the detector. But even in 
the central region, —1.5 < yiab < 1, to which previous measurements of charmed-meson 
photoproduction at HERA were conhned [0 , fragmentation production may be enhanced 
with respect to the fusion mechanism by selecting events with high photon energy. This 
is demonstrated in Fig. H, where the W dependence of (P'a/dW dpx is shown for pt = 5 
and 10 GeV. While the W distribution of the direct-photon contribution has a shape 
similar to that of the fusion contribution, with a maximum at low to intermediate W, the 
resolved-photon contribution tends to be monotonically increasing with increasing W. As 
a consequence, even at low pt, fragmentation production becomes more important than 
the fusion mechanism beyond some point in W. At px = 5, the cross-over takes place at 
about W = 125 GeV, which corresponds to x = 0.17 in Eq. ([T7|). 

Guided by this observation, we repeat the analysis of Fig. |l|(c) integrating over 150 GeV 
< W < 280 GeV, instead of 30 GeV < W < 150 GeV, and plot the outcome in 
Fig. As expected, this leads to a signihcant increase of the fragmentation contribution, 
which mainly originates from the resolved-photon part. On the other hand, the fusion 
contribution is only moderately increased at Pt ^ 4 GeV, while it is even decreased at 
lower values of pt- In conclusion, over the entire pt range, the fragmentation contribution 
is now more than twice as large as the fusion contribution; at GeV, the fusion 

contribution falls short of the fragmentation contribution by more than one order of 
magnitude. 

At HERA, it should also be possible to detect J/'ip mesons in the very forward di¬ 
rection, almost up to yiab = 3.5 [^. This would offer the opportunity to render the 
fragmentation to fusion ratio even more favourable. Thus, in Fig. |^, we investigate how 
the results for 150 GeV < W < 280 GeV in Fig. ^ are modihed if the yiab interval is 
shifted from —1.5 < yiab < 1 to 1 < yiab < 3.5. We observe that the fragmentation 
contribution is then enhanced, almost by a factor of two at large px- But, what is even 
more signihcant, the fusion contribution is reduced by one to two orders of magnitude. 
As a consequence, the fragmentation to fusion ratio is dramatically increased, to values in 
excess of 200 throughout the entire px range. We conclude, that the kinematical region of 
large W and large yiab is an ideal place to probe at HERA the anomalously large inclusive 
JI'll) yield observed at the Tevatron. 

It has become customary to present experimental data on inclusive J/'ip production in 
terms of the inelasticity variable f dehned in the Introduction, instead of using yiab- The 
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relation between 5 and yiab reads 


= 


2EpmT 


Vlab 


where Ep is the proton energy in the laboratory system, 
of interest corresponds to small values of 5. This is in line with Ref. 
kinematical cut 5 < 0.5 has been proposed to enhance the fragmentation contribution 
relative to the fusion contribution. We have recently assessed, in a NLO analysis a. 


(18) 

Thus, the kinematical range 
0 , where the 


the logistics of enhancing via appropriate acceptance cuts the colour-octet contribution 
to J/'ip production via fragmentation at HERA. 


4 Conclusions 

In addition to the well-known parton-fusion mechanism, J/ip mesons may also be produced 
via the fragmentation of hnal-state partons. The latter mechanism is expected to be the 
dominant source oi J/ip mesons at sufficiently large pt- In this paper, we studied the 
photoproduction oi J/ip mesons at large pt via fragmentation at HERA. 

Our calculation was performed at NLO in the QCD-improved parton model with 
Hf = A massless quark flavours, he., the charm quark was assumed to be an active flavour 
inside the proton and the resolved photon. We employed charm and gluon EE’s which 
were perturbatively calculated at a fragmentation scale of order Mj/^ in nonrelativistic 
QCD, and evolved them in 2 ; space to the characteristic scale of the considered process, of 
order px, using NLO Altarelli-Parisi splitting functions. At the same time, we adjusted 
the factorization scheme of the collinear singularities connected with hnal-state charm 
quarks so as to match the corresponding calculation with massive charm quarks. In this 
way, large logarithms of the type ln(py/Mjy^), which would be present in the massive 
calculation and render it unreliable at large px, are removed from the hard-scattering cross 
sections and properly resummed. We included prompt J/ip fragmentation production as 
well as Xcj fragmentation production followed by Xcj J/ip + x decays (non-prompt J/ip 
fragmentation production). In each channel, we took into account the contributions due 
to colour-singlet and colour-octet cc states. Apart from direct photoproduction, which 
had been considered to LO in Ref. [^, we also studied resolved photoproduction, both 
at NLO. 

While the nonperturbative matrix elements multiplying the colour-singlet EE’s may 
be reliably extracted from the decay properties of the charmonium states, the colour-octet 
matrix elements are adjusted to account for the enormous excess of recent Tevatron data 
on large-pr charmonium production over the theoretical prediction within the CSM. This 
procedure implicitly assumes, on the theoretical side, that there are no other produc¬ 
tion mechanisms that have been overlooked so far and, on the experimental side, that all 
backgrounds are fully under control. Therefore, the colour-octet matrix elements must 
be considered far less rigorously determined than the colour-singlet matrix elements. Ob¬ 
viously, it is of prime importance to probe the colour-octet matrix elements in other 
experiments as well. The signature for colour-octet J/ip production in e''“e“ annihilation 
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at the Z-boson resonance has recently been investigated in Ref. p6|,|3^. Here, we con 


tinned this research programme by assessing the prospects of the HERA experiments to 
substantiate the Endings at the Tevatron. 

Our main observations may be summarized as follows. As for direct J/^jj photoproduc¬ 
tion, fragmentation production starts to have a higher yield than the fusion mechanism of 
the GSM at px ^ 10 GeV [see Figs, ^(a) and @(a)]. This is in line with the main conclusion 
of the LO analysis of Ref. JT^, since the NLO corrections are relatively modest in this case 
[see Fig. 0(a)] . However, in the range px^lO GeV, which will be accessed at HERA in 
the near future, resolved J/V’ photoproduction via fragmentation, which has not yet been 
considered in the literature, is considerably more significant than direct photoproduction 
[see Figs. |I](b) and @(b)], by a factor of 18 (2.3) at = 3 GeV (10 GeV). In fact, the 
resolved-photon fragmentation cross section already exceeds the fusion cross section at 
Pt^5 GeV and is comparable to it at smaller values of px- The bulk of this cross section 
originates in the subprocesses + qp ^ g + X and q-y + gp ^ g + X followed by the 
fragmentation chain g —>■ cc[8, —*• J/'ip. The dominance of resolved photoproduction 

is partly due to a large QGD K factor, in excess of 2 [see Fig. 0(b)]. These results refer 
to the standard ranges —1.5 < yiab < 1 and 30 GeV< W < 150 GeV. 

We suggested two ways to further enhance the (resolved-photon) fragmentation pro¬ 
duction relative to the fusion mechanism. One is to concentrate on the forward direc¬ 
tion, e.g.^ 1 < yiab < 3.5, where the resolved-photoproduction cross section peaks [see 
Fig. 0(b)]. The other one is to increase the GM energy available for the hard yp scatter¬ 
ing, e.g., by choosing 150 GeV< W < 280 GeV, as the resolved-photoproduction cross 
section increases monotonically with W, while the fusion cross section takes its maximum 
at relatively low W values (see Fig. 0). The usefulness of increasing W and yiab is nicely 
illustrated in Figs. 0 and 0, respectively. Both options should be feasible at HERA, and 
we look forward to the exciting scrutiny of the Tevatron J/'ijj colour-octet phenomenon 
by the HERA experiments in the near future. 
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A Appendix: Coefficients of the perturbative frag¬ 
mentation functions 

For the reader’s convenience, we collect here the colour-singlet and colour-octet coefficient 
functions, at the starting scale po) for fhe gluon and charm-quark FF’s into J/il) and Xcj 
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mesons, with J = 0,1,2. The g —> cc[l, ^S*!] coefficient function has the two-dimensional 
integral representation 



where 


fo{r, y) = ^^(1 -I- r)(3 -I- 12r -|- 13r^) — 16r^(l -|- r)(l -|- 3r)|/ 

— 2r(3 — 9r — 21r^ -|- 7r^)y‘^ + 8r(4 -|- 3r -|- 3r‘^)y^ — 4r(9 — 3r — 4r^)j/'^ 

- 16(1 + 3r + 3r^)y^ + 8(6 + 7r)y^ - 32y\ 

/i(r, y) = —2r(l -|- 5r -|- 19r^ -|- 7r^)y + 96r^(l -|- r)y‘^ + 8(1 — 5r — 22r^ — 2r^)y^ 

+ 16r(7 + 3r)y^ - 8(5 + 7r)y^ + 32/, 
y) = ''"(1 + 5r -|- 19r^ -|- 7r^) — 48r^(l -|- r)y — 4(1 — 5r — 22r^ — 2r^)y‘^ 

- 8r(7 + 3r)/ + 4(5 + 7r)/ - 16/, 

9o{f'i y) = ■'"^(1 — 'r)(3 -|- 24r -|- 13r^) — 4r^(7 — 3r — 12r'^)y — 2r^(17 -|- 22r — 7r‘^)y‘^ 

+ 4r^(13 -|- 5r — 6r^)y^ — 8r(l -|- 2r -|- 5r^ -|- 2r^)y^ — 8r(3 — Hr — 6r‘^)y^ 

+ 8(l-2r-5r2)/, 

y) = —2r^(l -I- r)(l — r)(l -|- 7r)y + 8r^(l -|- 3r)(l — 4r)?/^ 

-|- 4r(l -|- lOr -|- 57r^ -|- 4r^)/ — 8r(l -|- 29r -|- 6r^)/ — 8(1 — 8r — 5r^)y^, 

92 ir, y) = r‘^{l + r)(l — r)(l -|- 7r) — 4r^(l -|- 3r)(l — Ar)y 

— 2r(l 4- lOr -|- 57r^ -|- Ar^)y'^ -|- 4r(l -|- 29r -|- 6r^)/ -f 4(1 — 8r — 5r^)/. (20) 

All other coefficient functions are available in closed form 





85 13 9 9 , ,, , 

- 1 + —z - -z^ + 7(5 - 3z) ln(l - z) 

8 4 4 ^ ^ ^ 






(l92 + 384z -t- 528^2 - 1376^^ + 1060z^ - 376z^ + 59z®) , 
(96 - 288z + 496;^^ - 408z^ + 202/ - 54/ + 7/) , 


2) + (/;)^ -l + ^ z - z ^ + ^( 2 - z ) 1,1(1 - z ) 


4/—T (le _ 32z + 12z? - 32z“ + bz *), 

243(2 - /6 V / 


(I 6 - 32^ + 72^2 - 32/ -f- 5z^) , 
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(z) = ^ f48 - 192z + 480/ - 668z^ + 54l/ - 184/ + 23/) . (21) 

3645(2 — z)^ ^ ^ 

As usual, the plus distribution 1/(1 — 2 ;)+is defined by dz f{z)/{1 — z)+ = Jq dz [f{z) — 

/(!)]/(! — z) for any regular function f{z). 

For the numerical solution of the /i^-evolution equations in z space, it is useful to 
approximate 


1 -z 


ln(l — z) = 


0 , if 2 ; < 1 — e. 


5(1-/ = 1 


-, if 1 — e < < 1 , 

e 

1 

, if 2 : < 1 — e. 


1 -z 


-r 1 

-, if 1 — e < z < 1 , 

e 

ln(l — z), if z < 1 — €, 

hie — 1 , if 1 — e < z < 1 , 


( 22 ) 


with appropriate 0 < e -C 1. We have verified, for the perturbative FF’s of Ref. |27 


m 


LO and NLO, that the /i^ evolution in space implemented using Eq. (^^ agrees very 
well with the /i^ evolution in Mellin space . 
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FIGURE CAPTIONS 

Figure 1: (a) Direct-photon, (b) resolved-photon, and (c) total contributions to the cross 

section da/dpT of inclusive J/ip photoproduction via fragmentation at HERA, integrated 
over 30 GeV < W < 150 GeV and —1.5 < yiab < 1- The NLO contributions due to 
prompt singlet (PI), prompt octet (P8), non-prompt singlet plus octet fragmentation 
(N1-1-8), and their sum (frag.) are compared with the LO GSM contribution (fusion). 

Figure 2: (a) Direct-photon, (b) resolved-photon, and (c) total contributions to the cross 

section da/dyiabPT of inclusive J/ip photoproduction via fragmentation at HERA, inte¬ 
grated over 30 GeV <W < 150 GeV, at = 5 and 10 GeV. The NLO contributions due 
to prompt singlet (PI), prompt octet (P8), non-prompt singlet plus octet fragmentation 
(N1-1-8), and their sum (frag.) are compared with the LO GSM contribution (fusion). 

Eigure 3: NLO to LO ratios for the (a) direct-photon, (b) resolved-photon, and (c) 
total contributions to the cross section da/dpT of inclusive J/ip photoproduction via 
fragmentation at HERA, integrated over 30 GeV <W< 150 GeV and —1.5 < yiab < 1- 
The ratios are taken for prompt singlet (PI), prompt octet (P8), non-prompt singlet plus 
octet fragmentation (Nl-l-8), and their sum (frag.). 

Eigure 4: Gross section d?a/dW dpr of inclusive J/ip photoproduction via fragmentation 
at HERA, integrated over —1.5 < yiab < 1, at pr = 5 and 10 GeV. The NLO fragmentation 
contributions due direct photons, resolved photons, and their sum are compared with the 
LO GSM contribution. 

Eigure 5: NLO fragmentation and LO GSM contributions to the cross section da/dpT 
of inclusive J/'^ photoproduction at HERA, integrated over —1.5 < yiah < 1 and two 
different W intervals, namely, 30 GeV < lU < 150 GeV and 150 GeV < IT < 280 GeV. 

Figure 6: NLO fragmentation and LO GSM contributions to the cross section da/dpx of 
inclusive J/ip photoproduction at HERA, integrated over 150 GeV < IT < 280 GeV and 
two different yiab intervals, namely, —1.5 < yiab < 1 and 1 < yiab < 3.5. 
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